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Abstract. We propose analytical lower and upper estimates on the excitation threshold for breathers (in 
the form of spatially localized and time periodic solutions) in DNLS lattices with power nonlinearity. The 
estimation depending explicitly on the lattice parameters, is derived by a combination of a comparison argument 
on appropriate lower bounds depending on the frequency of each solution with a simple and justified heuristic 
argument. The numerical studies verify that the analytical estimates can be of particular usefulness, as a simple 
analytical detection of the activation energy for breathers in DNLS lattices. 
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1. Introduction 

A great deal of attention has been paid to the study of localization phenomena in nonlinear discrete systems 
in recent years, interest which has been summarized in a number of recent reviews This growth has been 
motivated not only by its intrinsic theoretical interest, but also by numerous applications in areas as the 
nonlinear optics of waveguide arrays [2] , Bose-Einstein condensates [3] , micro-mechanical models of cantilever 
arrays [4] , or some models of the complex dynamics of the DNA [5] . 

In this framework, perhaps one of the the most prototypical model is the so-called, discrete nonlinear 
Schrodinger equation (DNLS) [HEI- DNLS may arise as a direct model, as a tight binding approximation, or 
even as an envelope wave expansion and, it could be possible to say that the DNLS is one of the most ubiquitous 
models in the nonlinear physics of dispersive, discrete systems [8]. 

Our aim in the present paper is to determine analytical lower and upper useful bounds for the formation of 
spatially localized and time periodic modes in focusing DNLS lattices, called discrete breathers or also DNLS 
solitons, with power nonlinearity [51 111)]. 

Flach, Kladko and MacKay [11], addressed the existence of energy thresholds for the formation of discrete 
breathers in in one-, two and three-dimensional lattices. They defined as energy thresholds, the positive lower 
energy bounds possessed by discrete breather families (DB). Their numerical findings and heuristic arguments 
considered a generic class of Hamiltonian systems and showed that the energy of a DB family has a positive 
lower bound for lattice dimension N greater than or equal to some critical dimension N c , whereas for N < N c 
the energy goes to zero as the amplitude goes to zero. 

For the focusing DNLS equation in the infinite lattice !* N , 

i^ n + e(A^)„ + A„|^„| 2ff -0, A>0, a>0, (1.1) 

where n = (711,77-2, . . . , njv) € 7L N , the hypothesis suggested by Flach Kladko and MacKay was resolved by 
Weinstein [12 . In (Ad^) n stands for the iV-dimensional discrete Laplacian 

(A d ?P) neZ N = V>m - 2A^„. (1.2) 

m£Af„ 

Here Af n denotes the set of 27V nearest neighbors of the point in Z N with label n. The parameter e > is a 
discretization parameter e ~ h~~ 2 with h being the lattice spacing and A > is the parameter of anharmonicity. 
The hypothesis of [llj was resolved in [12j for breathers in the ansatz of time-periodic solutions 

V» n (t) = e if2 Vn, ^>0, (1.3) 

spatially localized in the sense 

\ip n \ — * 0, as \n\ — > oo, 

(here \n\ — maxi<i<jv |Jii| for n = (ni,ri2, • ■ ■ , tin) & Z N ). 

Solutions (|1.3p of (jl.ip satisfy the infinite system of algebraic equations 

-e(A d 0)„ + ^„-A|0„| 2 ^„ = 0, n£Z N . (1.4) 

We can associate a power to any solution of the form (|1.3j) , defined as 

m = E i^i 2 - ( L5 ) 

The power (|1.5p together with the Hamiltonian 

H[4>] = e(-A d 0, 0) 2 - -L- V |0„| 2ff+2 , (1.6) 

(7+1 L — ' 

are the fundamental conserved quantities for (jl.ip . 

For the proof of [12] on the existence of the excitation threshold, a discrete version of a Sobolev-Gagliardo- 
Nirenberg inequality is crucial. This discrete version reads as 

E l^| 2CT+2 <^f E l^l 2 ) (~ A ^A)2, <7>|, (1.7) 

If C* is the infimum over all such constants for which inequality (|1.7p holds, then the excitation threshold 
^thresh is defined by [HI pg. 680, Eqn. (4.2)] 

(a + l)e(^thrc S h) _CT =C„ (1-8) 
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and the optimal constant C, has the variational characterization 
1 _ M (E„ e z« I<^IT(-A^) 2 



Then Weistein's result on the excitation threshold reads as follows: if 72. > T^thrcsh then I-ji < 0, and a ground 
state breather exists, that is, minimizer of the variational problem 

l n = inf{H[<j>] : V[</>] = H}. (1.9) 

On the other hand, if if 72. < 72-thresh then T-R = 0, and there is no ground state minimizer of (|1.9|) . In the light 
of the results of Weinstein, the critical dimension predicted by Flach, Kladko and MacKay is defined for the 
DNLS (dU) as 

N c = -. (1.10) 
a 

In this paper we propose analytical lower and upper estimates on the excitation threshold T^thrcsh, which 
are depending explicitly on the lattice parameters. The derivation of these estimates in section 2, can be briefly 
described as follows: for the lower estimate, we derive first by a fixed point argument, a lower bound on the 
power of the breather solution satisfied for any SI > 0. The role of such local bounds (through their dependence 
on the frequency Q) as thresholds on the existence of breather solutions has been analyzed in detail and tested 
numerically in [131 114] . Then this is compared with a second local lower bound involving this time the unknown 
value of 72.thrcsh- Although the lower bound for T^thrcsh derived as above, depends on an unspecified positive 
integer, its appropriate value can be easily determined by a simple and justified heuristic argument, explained 
in detail in subsection 2.1. The derivation of the upper bound comes out by simply examining the interpolation 
inequality (|1.7|) in comparison with the standard embedding inequality between the ^-sequence spaces. 

The numerical studies performed in subsection 2.1, justify that the estimates for 72.tincsh can be useful 
(on the account of their explicit dependence on the lattice parameters and the simplicity of the formulas), 
in "trapping" the exact value of 72-throsh for the cases of nonlinearity exponent a and dimension N which 
are of primary physical interest. This "trapping" is of particular interest in applications since the analytical 
estimation of the excitation threshold can be used for a simple calculation of the activation energy needed for 
the experimental detection of discrete breathers [11] . It is important to recall that the excitation threshold 
appears in the formal continuum limit e — > oo only in the case a — 2/N |12j . 

2. Analytical lower and upper bounds for 72.thresh- 

Our arguments on the determination of simple analytical bounds for the excitation threshold 72. t hrcsh will be 
based on some technical lemmas involving lower bounds for the power of solutions (|1.3|) for all Q > 0. 

Lemma 2.1 The power of a nontrivial breather solution of satisfies the lower bound 

i 

O 

_4eA7V(cr + 1) 

Proof: Multiplying (11.41) in the ^ 2 -scalar product we infer that </> satisfies the energy equation 

e(-A d cf>, 4>) 2 + nJ2 \^\ 2 =^J2 \^\ 2a+2 , for all fi > 0. (2.12) 



^-min.l ■ ^-thresh 



<n[4>] for all n>0. (2.11) 



Now inserting the inequality (|1.7p in the right-hand side of (|2.12p . and noting that 

HW,0) a <4iV Y l^l 2 ' ( 2 - 13 ) 

nei N 



we deduce that 



e(-A d 0,0) 2 + fl I0«| 2 < AC, ( \<t>nA (-A d 

neZ N \n£Z N / 
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ft E I^l 2 < AkNC * ( E l^l 2 ) 



tr+1 

(2.14) 



By substitution of (| 1 . T[) into (|2. 14[> . we derive the lower bound (|2.1ip . o 

Lemma 2.2 Let k € M + , k > i, arbitrary. Then every non-trivial breather solution fl.3\) of 11. 1)) has power 
satisfying 

V2k- 1 ft 



< TZ[4>] for all ft > 0. (2.15) 



k A(2a + 1) 

Proof: We recall and modify appropriately the fixed point argument of [13l[T5]. We consider the operator 

- eA d + ft : -> £ 2 . (2.16) 

being linear and continuous. It also satisfies the assumptions of Lax-Milgram Theorem [TBI Theorem 18. E, pg. 
68]: Note that 

e(-A d «^) 2 + ft||«£||! >tl\<t>\\l forall^e^ 2 . (2.17) 
Then according to Lax-Milgram theorem, for given z 6 £ 2 , the linear operator equation 

-eA^ n + n^ n = A\z n \ 2<T z n , A>0, (2.18) 
has a unique solution <fi £ £ 2 , since 

w\z\ 2 °z\\i< e ^r +2 <iNr +2 - 

Hence we are allowed to define the map V : I 2 — > ^ 2 , by 'P(z) := where </> is the unique solution of the operator 
equation (|2.18[) . Clearly the map V is well defined. Let £,£ be in the closed ball 

B fl :={ze< 2 : \ \z\\p < R}, 

such that <j) = ^(C)j V' = "P(Q- The difference x := — ^ satisfies the equation 

- eA dX „ + ftXn = A(|C„rCn - |6»| 2 *£») (2.19) 

We recall that for any F € C(C, C) which takes the form F(z) — g{\C,\ 2 )Ci with g real and sufficiently smooth, 
the following relation holds 

F(C)-F(0= [ {({-0(9(r)+rg'(r)) + ((-l)<S> 2 g'(r)}dd, (2.20) 
Jo 

for any C, £ G C,where $ = 6>C + (1 - 0)£, 6> e (0, 1) and r = |$| 2 (see [T7| pg. 202]). Applying (|2T2Uj) for the 
case of F(C) = |C| 2<7 C one finds that 

ici 2ct c - = [ i(a+i)(c- om 2 ° +<?(<- o&m 2 °- 2 }d&. (2.21) 

Jo 

Assuming that £, £ s Br, and noting that ||$||2 < -R, we get from (|2.21j) the inequality 

E nc»rc« - I6,re„i 2 < (2a + 1) 2 E { f i $ «i 2c7 ic™ - 

<(2a + l) 2 E f ril^lla'Kn- Snide! 

< (2a + 1) 2 E ( I' R 2a \(n-Udo) 

= (2a + l) 2 i? 4 " E K«-a 2 - (2-22) 
Taking now the scalar product of (|2.19|) with x m ^ 2 an d using (|2.22|) . we have 

<A||xll2|||CI 2ff C-|£l 2<T £l| 2 

<A(2a + l) 2 J R 2CT ||xll2||C-CI| 2 . (2.23) 



e (-A d x,x)2 + ft||xll^<A||x|| 2 |||C| 2CT C-|e| 2<T eil2 

2 d 2ct| 



Lower and upper estimates on the excitation threshold for breathers in DNLS lattices 
Applying next Young's inequality 



6 v 1 
ab < -aP H j- 

p qei/P 
withp = o = 2, a= ||x|| 2 , & = ||C-£I|2 and 

e=-, k>1/2, 
in (HH]), we get that 



o 9 , for any e > 0, l/p+l/q=l, 



(2k- l)fi 
2^ ' 



X||l< ^A 2 (2 ( T+l) 2 i? 4 -||C-C||i. 



(2.24) 



From (|2.24|) . we conclude with 



IWll = ITO-n0ll2< 



-A 2 (2a + l) 2 i? 4CT ||C-?ll^ 



ft 2 (2fc- 1)" ' iy " s " 2 ' 

Since 'P(O) = 0, from inequality (|2.24p we derive that the map V : Br — > is a Lipschitz map with 
Lipschitz constant 



L = 



=A(2ct + l)i? 2 



f7V2fc - 1 

The map V will be a contraction and will have a unique fixed point if L < 1. This unique fixed point will 
be the trivial one, since V(0) = 0. Hence, for 



R 2 < 



V2K- 1 



A(2o-+l) 



the only breather solution is the trivial. Therefore, a non-trivial breather solution (jl.3[) should have power 

> ^mi„,2- O 

Keeping the positive constant n > 1/2 undetermined is crucial to derive a K-dependent lower bound for 
^•thresh- However the numerical investigations, will reveal that for practical purposes the constant k can be 
easily determined. We start by proving the following 

Proposition 2.1 Let a > 2/N. There exist K cr a > 1/2 such that 



V2Kcrit - 1 4Ae((7+l) 



<-Rthrcsh < [4eN(<r + !)]'• 



(2.25) 



Kcrit 2a +1 

Proof: It follows from Lemma \2. II that lirn K ^oo TZ m i n ^(K) — 0, therefore, we can make 7?. m i n| 2(/c) as small as 
we please by taking k large enough. Thus, there exists a K cr it such that 7?.min,2(K C rit) < 7?-min,i, implying the 
left-hand side of ()2.25|) . For the right-hand side of (j2.25[> . we recall first that 



E W ^ E \^\" , for all 1 < g < p < oo. 



(2.26) 



ne'i 



Applying the inequality (|2 . 26[) for p = 2a + 2 and <; = 2we get that 



E \<t>« 



|2ct+2 



< 



E 



for all cr > 0, (f>e £ 2 . 



From (|2~271) we have 



sup 



2cr+2 



ee 2 (E„ 6Z «l0nl 2 ) 



—rr < 1, for all a > 0. 

r-j-l — — 



Since the above inequality holds for all a > and € £ 2 , by setting = e J for arbitrary j £ N, any element 
of the orthonormal basis of £ 2 , in (|2.26j) , we get that 



sup 

(pee 



2 E„ e Z«l^| 2CT + 2 



1, /or all a>0. 



(2.27) 
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On the other hand, it follows from (jl.7p and (|2. 13|) that 



Vn6Z« / 



(7 > 2/N, 



which implies that 

(SngZ^ l^ l 2 ) < 4ArC * for aU a > 2/Ar g £ 2 



v cr+l 



(2.28) 



Then a comparison of (12.27| which holds for all a > 0, with (|2.28|) implies that 

1 < WC* - 4e/V(<7 + l)iZ^, (2.29) 

from which we conclude the right-hand side of (|2.25[) . Even in a simpler way, one can set in (|2.28[) any element 
of the orthonormal basis of £ 2 to derive (|2.29[) . o 

3. Numerical study 



Since the upper bound on (|2.25[) is exact, the estimates would have a full strength in applications, if the 
undetermined constant K cr ;t could be easily determined, at least by a simple heuristic argument. For such a 
simple heuristic determination of the constant /c C ritj it looks natural to restrict to the case k G Z + , k > 1. Then, 
the simplicity of the formula (|2.25[) suggests that the appropriate value of n can be determined by considering 
successive choices of k. Setting 



and rewriting (|1.8[) as 



V2Kcrit - 1 4iVe(<r+l) 



^crit 



2(7+1 



^•thresh 



(a + l)e 



we observe that 



lim 1Zr, 



lim ^thresh = I- 



(3.30) 



independently of the choice of k, e, N . This behavior completely justifies that even the first choice K cr j t = 1 is 
valid for "sufficiently large" a. The fist numerical study whose results are demonstrated in Figure [TJ examines 
the range of a > on which this simplest choice K cr it = 1 is valid, i.e. the validity of the formula 



2(7+1 



< -Rthrcsh < [4dV(tr+l)p 



(3.31) 



The green dashed line represents the theoretical upper estimate lZ u b := [AeN(a + l)p, the blue full line 
corresponds to the numerical -Rthresh as a function of a > 2/iV and the red dashed line represents the theoretical 
lower estimate IZib- The first numerical study, not only reveals that the formula (|3.31[) is valid for the case 
N = 1,2 but also of very good accuracy for N — 2 and excellent for N = 3 for a > 1 with a discrepancy 
regarding the prediction of the lower bound TZib appearing in the interval a £ (2/3,1). In the light of the 
behavior i3.30\) . the choice K cr j t = 1 is satisfied for all a > 1. Motivated by the recent work of J. Dorignac, 
J. Zhou and D.K. Campbell [H] which considers integer values of a > 2/N (represented by the black dots in 
the figures) it seems fair to state that the prediction of (|3 . 3 1 1) is of particular usefulness for such nonlinearity 
exponents and lattice dimensions which are of main physical interest. Seeking for the value of K cr it which 
would remove the small discrepancy of (|3.31j) for N = 3 and real values of a > 2/N, our numerical findings in 
Figure [2] verified that in the choice K cr ;t = 2 this discrepancy is reduced to the interval a € (2/3,0.72) and it 
is completely removed for the choice of K cr ;t = 3, as it is shown in Figure [3J A summary of our findings for 
the cases N = 1,2,3, suggests to restate Proposition 12.11 taking into account the dependence of K crit on the 
dimension of the lattice: Letting n € Z + and N being fixed, we observe that since liniK^oo TZ m in,2{^) — 0,we 
can always find K CI it(N) > TV such that 



y/2 Kcrit (N) - 1 4iVe((7 + l) 



Kcrit (N) 



2(7+1 



< Rthresh < [4e/V((7 + 1)] ' , for all 1 < N < K crit (/V).(3.32) 
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Figure 1. Numerical values for 7^ t hresh as a function of a > 2/N against its lower and upper estimation l|2.25|l 
for re crit = 1 (formula (OT}). (a) N = 1, a > 2, (b) N = 2, a > 1, (c) N = 3, cr > 2/3. In all cases e = 1. Green 
dashed line corresponds to the upper estimate, blue full line to the numerical 7£ t hresh an <i re d dashed line to 
the lower estimate The inset in (c) magnifies the discrepancy observed for the prediction of the lower estimate 
of l|3,31| l in the interval a £ (2/3, 1). Black dots correspond to integer values of the nonlinearity exponent a. 



With the rigorously valid estimates (|3 . 32(1 at hand, the numerical study for the cases N = 1,2,3 suggest that 
when N = 3 it is justified to consider K C rit(^V) = 3 and that 



y/5 Weja + 1) 
~3 2cr+l 



< -Rthrosh < [4e7V((j + 1)] * , for all 1 < N < 3. 



(3.33) 



Actually the numerical study and especially the collection of Figures [T] (a), (b) and [3] (c) justify the 
validity of the formula 



ViV^T 4Ne(a+l) 



N 



2(7+1 



< ^thresh < [4eN(o- + l)p , for all 1 < N < 3, 



(3.34) 



which is of valuable accuracy for N = 2,3. The estimates (I3.33|) . (|3.34p have the advantage of removing the 
small discrepancy of (|3 . 3 1 [) observed in the case N = 3, for real a > 2/N. However we believe that all the above 
formulas derived by a simple heuristic implementation of Proposition ^. 1[ serve as a very satisfactory analytical 
estimation of the excitation threshold in the cases of cr, N which are of physical significance. 
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Figure 2. Numerical values for 7^ t hresh as a function of a > 2/7V against its lower and upper estimation l|2.25|l 
for K cr i t = 2. The inset in (c) magnifies the discrepancy observed for the prediction of the lower estimate of 
l|3.32ll in the interval a £ (2/3,0.72) which is reduced in comparison with Figure [T] (c) . 



4. Conclusions 

In this work, we have determined analytical upper and lower estimates on the excitation threshold for breathers 
in TV-dimensional DNLS lattices. Numerical calculations show that, in cases studied, the theoretical bound is 
close to the true threshold providing useful analytical expressions to determine analytical energy activation of 
breathers in these systems. On the other hand, extensions of previous results to more general situations, as 
DNLS systems with impurities, are currently under investigation and will be reported in future publications. 
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